Quantum rings with time dependent spin-orbit coupling: 
Rabi oscillations, spintronic Schrodinger-cat states, and conductance properties 
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The strength of the (Rashba-type) spin-orbit coupling in mesoscopic semiconductor rings can 
be tuned with external gate voltages. Here we consider the case of a periodically changing spin- 
orbit interaction strength as induced by sinusoidal voltages. In a closed one dimensional quantum 
ring with weak spin-orbit coupling, Rabi oscillations are shown to appear. We find that the time 
evolution of initially localized wave packets exhibits a series of collapse and revival phenomena. 
Partial revivals - that are typical in nonlinear systems - are shown to correspond to superpositions 
of states localized at different spatial positions along the ring. These "spintronic Schrodinger-cat 
sates" appear periodically, and similarly to their counterparts in other physical systems, they are 
found to be sensitive to environment induced disturbances. The time dependent spin transport 
problem, when leads are attached to the ring, is also solved. We show that the "sideband currents" 
induced by the oscillating spin-orbit interaction strength can become the dominant output channel, 
even in the presence of moderate thermal fluctuations and random scattering events. 

PACS numbers: 85.35.Ds, 03.65.-w 
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I. INTRODUCTION 

The observation of fundamental flux- and spin- 
dependent quantum interference phenomena that 
can appear in quantum rings made of semiconduct- 
ing materials 1 exhibiting Rashba-type^ spin-orbit 
interaction 3 -^* 5 - (SOI) have motivated many studies in 
the past few years. Additionally, as the strength of 
the SOI that determines the spin sensitive behavior 
can be tuned with external gate voltages^, quantum 
rings or systems of them&i&^iiii also have possible 
spintronic applications. Because of the conceptual 
importance and the possible applications, closed single 
quantum rings (without attached leads) ) 11 ' 12 ! 13 ! 14 ! 15 
as well as two- or three-terminal ones were 
invectigate d 15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23,24,25,26 - 27,28,29 - 30,31 - 32 - 33,34 
extensively. 

In the current paper we explore effects which are ex- 
pected to appear when the SOI in a single quantum ring 
is time dependent. The strength of this interaction is as- 
sumed to be a sinusoidally oscillating function of time. 
Studies of transport related problems with oscillating 
SOI have been initiated in Ref. [35| for a ring, and in 
Ref. [36| for a ring-dot system, mainly in the context 
of spin currents. Here we focus on different aspects by 
starting the analysis with the case of closed rings and de- 
termine the energy levels and eigenstates of the relevant 
Hamiltonia n 37 ! 38 . Using these analytic results we can cal- 
culate the dynamics for arbitrary initial conditions. The 
time evolution of states with well-defined total angular 
momentum component j in the direction perpendicular 
to the plane of the ring is shown to be analogous to the 
classical Rabi flopping. That is, the spin components os- 
cillate with the Rabi frequency, which in our case depends 



on the eigenvalue j, and the oscillation itself is relatively 
stable against disturbances caused by random scatterers. 
When the initial state is a localized wave packet, then, 
theoretically, the time evolution is periodic, if the ratios 
of the relevant eigenfrequencies are rational. We show 
that appropriately chosen amplitude of the SOI oscilla- 
tions leads to quasi-periodic dynamics, even when the 
relevant frequencies are not exactly commensurable. 



The description of the conductance properties of the 
ring requires the solution of a quantum mechanical scat- 
tering problem with a time dependent Hamiltonia n 35 ' 36 . 
In this case energy is clearly not a constant of motion, 
the relevant continuity equation contains a source term 
describing the energy explicitly calculate this source term 
as well as the energy current density. These results pro- 
vide a clear physical picture from which we can solve the 
scattering problem. As an application, we show that the 
ring can shift the energy of the incoming plane waves 
(both up- and downwards) by the frequency of the SOI 
oscillations expressed in energy units. This effect means 
the emergence of "sideband" currents in the transmis- 
sion, which, in our case, can become stronger than the 
direct one. Clearly, harmonics of the driving frequency 
(which corresponds to the sideband currents in our case) 
appear naturally in driven nonlinear systems. However, 
let us emphasize that quantum rings are special in the 
sense that their geometry induces nontrivial interference 
effects that determine the energy dependent transmis- 
sion probabilities. This characteristics of the device is 
responsible for the increased relative weight of the side- 
band currents in the output. 
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II. SPIN OSCILLATIONS IN A CLOSED RING 
A. Model 

We consider a ring 3 -^ of radius a in the x — y plane and 
assume a time dependent electric field in the z direction 
controlling the strength of the spin-orbit interaction char- 
acterized by the parameter o£. The Hamiltonia n 37 i 38 in 
the presence of spin-orbit interaction for a charged par- 
ticle of effective mass m* is given by 



h = nn 



.d_ 

'dip 2VL ' 



An 2 



(i) 



where ip is the azimuthal angle of a point on the ring (sec 
Fig. [IJ and the radial spin operator is given by a r /2 = 
(a x cosip + (TySmip)/2. Additionally, Ml = h 2 /2m*a 2 
denotes the kinetic energy of the charged particle and 



uj(t) = a(t)/ha = 2Acos(i4), 



(2) 



with v being the circular frequency of the periodic exter- 
nal electric field. Nanoscale quantum rings, for which 
the Hamiltonian above is relevant, can be fabricated 
from e.g. InAlAs/InGaAs based heterostructures 4 ^ or 
HgTe/HgCdTe quantum wells 1 -. For a ring of radius 250 
nm made of InGaAs, the frequency corresponding to f2 
is around 10 GHz. 
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FIG. 1: The geometry of the device and the relevant spinor 
valued wave functions in the different domains. 

In a closed ring, when there are no leads attached 
to it, the spinor valued wave functions have to be pe- 
riodic in space. Thus we look for the solution of the 
time evolution induced by the Hamiltonian ([T]) in terms 
of states which are eigenvectors belonging to an integer 
eigenvalue n of the z component of the orbital angular 
momentum L z — —id Vl and also eigenvectors of the spin 
component S z — er 2 /2 (in units of h). These operators, 
however, do not commute separately with H, whereas 
in each fixed moment of time we have the commutator 
[H, L z + S z ] = 0, therefore an eigenvalue of J = L z + S z 
is a good quantum number—. A given eigenvalue j of 
J is, however, still doubly degenerate, the two space- 
dependent spinors 



lT,n> 







belong to the same: j = n + 1/2 half integer eigen- 
value of J. As a consequence - similarly to the case of 
constant uj - the eigenvalue equation of H expanded in 
the {|Tji) ■ II) ti + 1)} manifold reduces to separate 2x2 
matrix problems corresponding to given values of j. Re- 
stricting H to one such subspace we obtain the two- 
dimensional matrix 



Hi 



fn(j - 1/2) 2 



Mt) 

Q(j + 1/2) 2 



In order to solve the time dependent Schrodinger equa- 
tion with the periodic coupling of the form ([2]) we look 
for the solution in a subspace fixed by a given value of j 
in the form 

\m) - ^{ty^ it, j - iM+b&y-w + 1/2) . 

(3) 

This is useful when - as in the current paper - the focus is 
on relatively weak spin-orbit strengths, and one can use 
standard rotating wave approximation (RWA), i.e., terms 
oscillating rapidly, with frequency 2v in the off-diagonal 
terms in the Hamiltonian can be neglected. Note that 
this widely used approximation is essentially equivalent 
here to the first order corrections (with respect to v) in 
an appropriate Floquet scattering matrix descriptio n 4 ^ 43 
of the problem. Using RWA, and introducing the dimen- 
sionless time variable r = fit, the system of equations 
determining the evolution has the form: 



H 4 



with the dimensionless operator 



3 bi 



U 



(4) 



(5) 



The solution of the system ([J| amounts to solve the 
eigenvalue equation of Hj yielding 



Ef = (f + ±)±y(2j-i>) 2 + l \\r. ((>) 



where v = v/Q.,A = A/Q. The corresponding eigen- 
spinors are given by 



(7) 



where Uj — cos9j/2, Vj — sin 6j/2, tan^- = ~ 2 J 4 ^. and 



u.) 



(2j 



-n2 
V) 



4A 2 j 2 + 2j 



2Aj 



(8) 



Note that for a given j, the spin operator S(6j,ip) — 
S x sin 0j cos tp + S y sin 6j sin tp + S z cos 9j (the spin com- 
ponent in the direction given by 9j and tp) also commutes 
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with H, therefore in the case of a given j the expecta- 
tion value of the spin direction varies along the ring in 
accordance to the ip dependence of the operator above. 

The states {|^>,i = ... -3/2,-1/2, 1/2,3/2...} 
obviously form a basis in the space of periodic spinor 
valued wave functions on the ring. Therefore, the time 
evolution (in the interaction picture introduced above) of 
any initial state 



(m 



(9) 



can be obtained in a straightforward way: 

iv>w> = E ( e " £j+T K f > (i^t i ( 10 ) 

3 

+ e-^|^>((^-|V(Q)»), 

where the inner products ((.|.)) that provide the expan- 
sion coefficients are the following: 



1 



«^ + I ^(0)» = u*^ j e-^-^f&W 
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v*^ j e- i « +1 ™*g(<p)d<p.(ll) 



B. State evolution of electrons with definite 
angular momentum: Rabi oscillations, collapse and 
revival 

States with well defined z component of the total an- 
gular momentum are linear combinations of the spinors 
\ipf) and \ipj) with the same value of j. The time evo- 
lution of these states can be of interest because conser- 
vation of the angular momentum may provide a method 
for preparing them. Absorption of circularly polarized 
photons e.g. can excite these states. (Interaction of short 
light pulses with the electrons confined in a ring has been 
discussed in Refs. [44M5I ].) 

As it is knowr4£, if a two-level system is getting excited 
with a resonant external field, the interaction produces 
Rabi oscillations of the level populations with a frequency 
proportional to the amplitude of the external field. If 
there is a detuning between the external field and the 
level splitting, the Rabi frequency is modified and the 
oscillations become less pronounced. The dimensionless 
Rabi frequency corresponding to the Hamiltonian Hj (see 
Eq. ((S])), is given by 



(2j -vf +4A 2 f. 



(12) 



In our case the (dimensionless) level splitting AE = 2j, 
the detuning 2j — z>, as well as the effective coupling Aj 
depend on j, that is, on the z component of the total 
angular momentum. This implies that the resonance fre- 
quency of the external field is also different for different 



values of j, namely v r = v r fl = 2jtt. When this resonance 
condition is met, the weights of the spin up and spin 
down components in the eigenspinors given by Eq. ([7]) 
are equal, tan Oj/2 = 1 (in this case independently from 
j.) Assuming that the initial state is spin-polarized in the 
positive z direction, 







(13) 



complete Rabi oscillations appear in the resonant case. 
That is, the time dependent quantum mechanical expec- 
tation value of the z component of the spin oscillates 
between —1/2 and 1/2 (in units of K) at a given point of 
the ring: 



S z (r, ip) =S Z (0, ip) cos n R T 



1 



■ cos Qrt. 



(14) 



(The amplitude of these oscillations are smaller for non- 
resonant external fields.) In the time evolution of the S x 
and S y expectation values, oscillations with the frequency 
of the external field superimpose on the Rabi flopping: 



Sx(t, ip) 
S v (t, tp) 



1 



sin(riflr) sin(j>r — (p) 



(15) 



1 



sin(f2flr) cos(£r — tp). 



Fig. [U visualizes the spin precession for the resonant case. 
An arrow starting from a certain point of the ring points 
into the direction of the spin at that spatial positio n 12 ' 47 . 
The lengths of the arrows are proportional to the local 
electron density regardless of the spin direction: 



pfat) = (¥(p,t)|¥fo>,t)>, 



(16) 



with (.|.) denoting the inner product of spinors (without 
integrating over the spatial degrees of freedom.) Note 
that in the current case p - and consequently the lengths 
of the arrows - do not depend on the position, in this 
sense the oscillations affect only the spin degrees of free- 
dom. 

Let us note that if the initial state is not a pure quan- 
tum mechanical state, e.g., it is an equal weight inco- 
herent sum of the eigenstates (that is, the spin part of 
its density operator is proportional to unity), then no 
Rabi oscillations will be visible. Therefore in a usual, 
not specially prepared sample in equilibrium at low tem- 
peratures the question whether these oscillations can be 
observed depend on the position of the Fermi level: if it is 
between E^ and EJ for some value of j, then Rabi flop- 
ping can be present, while when all occupied states are 
"paired", i.e., the eigenstates corresponding to Ej and 
EJ are both occupied, their opposite precession cancels 
out. 



C. Time evolution of wave packets 

In contrast to the previous subsection, where we have 
treated states with well defined angular momentum, we 
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FIG. 2: Resonant Rabi oscillation in the quantum ring. The 
initial wave function is given by Eq. (|13[) with j = 3/2, while 
A = 0.1, leading to D.r = 0.3. See EPAPS Document No. 
[number will be inserted by publisher] for the related movie 
file. 



consider here the time evolution of a complementary sit- 
uation, when the initial state is a localized wave packet 
with a small uncertainty in the position variable ip. Start- 
ing with a narrow, spin-polarized initial state, one ex- 
pects it to diffuse so that the spin direction changes 
locally during the process. However, the ring geome- 
try implies that the "tail" and "head" of the spread- 
ing wave packet will interfere when they start to over- 
lap. This leads to rather complex dynamics with an ad- 
ditional characteristic feature that can be seen at this 
point. Namely, if we assume that there is a finite number 
of non-negligible coefficients in the expansion of the ini- 
tial state in terms of the eigenstates of the Hamiltonian, 
the discrete nature of the spectrum (which is again a con- 
sequence of the geometry) may cause the initial phases to 
be restored after a certain "revival time" . In other words, 



we expect periodic "collapse and revival" phenomena: 
The initially localized wave packet becomes delocalized 
along the ring, but later it reassembles again. (Note that 
the term "collapse" - similarly to the case of an atom 
interacting with a quantized field 4 ^ - means here merely 
the decay of some expectation values and completely un- 
related to the notion of measurement induced " collapse of 
the wave function" .) In principle, when the ratio of some 
important frequencies is irrational (i.e., it is not a fraction 
of two integers), the revival time is infinitely long. How- 
ever, revivals appear when all frequencies can be written 
essentially as integer multiples of a base frequency. This 
condition can be met by appropriately choosing the am- 
plitude A = A/Q of the SOI oscillations in Eq. |2]) so 
that 1 + A 2 = (m/k) 2 with integer k and m. Then 

E n ~{f + fc,m integer, (17) 

and the base frequency is Q/k, leading to revival times 
T r = 2kn/fl. The approximation above is valid if v is 
negligibly small compared to the relevant values of j. 
Clearly, this requirement cannot be met always, but as 
we will see, rapidly oscillating SOI strengths modify the 
above picture only in the sense that revivals become less 
pronounced. 

Fig. [3] shows the time evolution of the electron density 
given by Eq. (|16[) at ip = for a wave packet which is ini- 
tially polarized in the positive z direction and centered at 
ip = with Gaussian envelope (see Fig. |U) In the upper 
(lower) panel A 2 = 3 (21/4), Eq. (T7]) is satisfied with 
m = 2,k = 1 (m — 5,k — 2), thus the revival time T r 
is very close to 2tt/H (Att/H). The initial phase relations 
are restored periodically, although for larger values of v 
the amplitude of the consecutive revivals decay faster. 
Clearly, the requirement of commensurable frequencies 
is related to a rather mathematical point of view, it is 
hardly possible to exactly satisfy it in an actual exper- 
iment. Therefore we investigated less ideal parameters 
as well, and, according to the inset of Fig. [3J less pro- 
nounced revivals appear also when the parameters are 
not exactly the ideal ones. 

However, an additional, genuinely nonlinear effect can 
also be seen in Fig. [3J partial revivals at certain time in- 
stants T r /m, with m being an integer. The significance of 
these less pronounced peaks can be seen in Fig. [H where 
the dynamics of the spinor valued wave function is visu- 
alized in the same way as in Fig. [2j As we can see, the 
initially localized wave packet becomes first delocalized 
(collapse), then, at the partial revival times there is a 
superposition of wave packets localized at different posi- 
tions, and finally, at T r , we can see a single wave packet 
again (revival). Let us note that the emergence of the 
"spintronic Scrodinger cat" states is a typical nonlinear 
feature, similar phenomena appear e.g., in the case of 
a wave packet moving in a Morse potential 49 . Finally, 
let us note that so far we considered only the effect of 
the energy levels, but the oscillation of the external field 
induces an additional, overall rotation of the spin direc- 
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FIG. 3: The time evolution of the electron density at ip = 
for the initial state shown in Fig. [4] Th e para meters are A = 
\[%,v = 1.0 (upper panel) and A = ^/21/4, v = 0.1 (lower 
panel); these values satisfy Eq. (|17|l with m = 2, k = 1 and 
m = 5, k = 2, respectively. The difference in the periodicity 
can clearly be seen, and signatures of partial revivals are also 
present. The inset of the first panel shows the dependence of 
the height of the first revival peak (denoted by the arrow) as 
a function of A. 



tion around the z axis at a frequency of v. This is why 
although the spatial form of the initial wave packets are 
restored at T r , the spin direction is usually different from 
the initial one. 
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FIG. 4: Wave packet motion in a quantum ring (A = v = 
1.0.) Note the appearance of "spintronic Schrodinger-cat" 
states. See EPAPS Document No. [number will be inserted 
by publisher] for the related movie file. 



III. RINGS WITH ATTACHED LEADS: 
CONDUCTANCE PROPERTIES 

A. Local energy balance 

One of the most direct ways of gaining information 
of a semiconductor device is measuring its conductance. 
Clearly, the model discussed so far cannot predict prop- 
erly the result of such experiments, effects related to the 
leads that connect the device to the contacts have to be 
considered as well. When the Hamiltonian of the ring 
does not depend on time, this problem has already been 
solved^. In that case basically three requirements should 
be met to obtain the solution: i) Energy conservation, ii) 
continuity of the wave function, and iii) vanishing net 
spin current densities at the junctions. Calculations of 
the conductance are usually carried out at a fixed energy 
(practically at the Fermi energy, as at low temperatures 
electrons around this energy level determine conductance 
properties) , in which case energy conservation is trivially 
satisfied. Having solved the eigenvalue problem of the 
Hamiltonians in the different domains (ring sections and 
leads), points ii) and iii) above imply that the resulting 
wave functions have to be joined together. The spinor 
components have to be continuous, and the spin current^ 
that enters a junction, should also leave it: 



all leads 



(18) 



Let us note that besides the Griffith's boundary 
conditions^ that have been described above (and will 
be used throughout this paper) there are other physi- 
cally realistic and often used possibilities as well. The 
choice of the boundary conditions is in fact shown to be 
related to the reduction of a two dimensional problem to 
one dimension 51 . 

However, in our case, when the Hamiltonian itself con- 
tains explicit time dependence, energy conservation can- 
not simply be taken into account by solving the problem 
within a given eigensubspace of the Hamiltonian. In- 
stead, we use a continuity equation containing an explicit 
source term S 

-p B (v,t)=Vj E +S=—j E (<p,t)+S(ip,t), (19) 

for the local energy density 

p E fat)=RB(V{<p,t)\H(<p,t)*(<p,t)), (20) 

with the spinor inner product of Eq. (JTHJ). Note that, 
however, the Hamiltonian ([1]) is Hermitian only with re- 
spect to the inner product that involves spatial degrees 
of freedom as well (see Eq. (fTTj) ). that is why we take the 
real part of the expectation value in Eq. (|20[) . Let us note 
that the time dependent Hamiltonian ([I]) can be inter- 
preted to describe the circular motion of a nonrelativis- 
tic charged particle in the presence of the scalar poten- 
tial hio 1 4ef2 and a spin dependent vector potential. This 
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latter means an effective, spin-dependent electric field 52 
being proportional to the time derivative of ui(t), and can 
also be related to an effective electromotive forced that 
can induce spin currents. 

A calculation similar to that given in the Appendix of 
Ref. [H[ leads to: 

S(cp,t) = Re<*fo>,t)|J^%, (21) 



and 



3e{% <P) 



Re 



ifi((*| — HV) 

dip 



'dip 



(22) 

The continuity equation ([T9j) is a local relation, its phys- 
ical meaning is seen most clearly when it is integrated 
over a certain domain: the change of the energy inside 
the domain is a consequence of the energy currents that 
flow in/out through the boundaries, plus the source term 
related to the time dependence of the Hamiltonian. In 
our case this is due to the oscillating strength of the spin- 
orbit interaction, so basically the time dependent electric 
field between the gate electrodes provides this extra en- 
ergy. In the limit when the domain reduces to a single 
point, e.g., to one of the junctions, finite terms (e.g. the 
source term) that are integrated over this domain of zero 
measure disappear, thus Eq. (fT§|) reduces to 



(23) 



E 3e = 0, 



all leads 

that is, the net energy current density has to vanish. 



B. Solution with time dependent boundary 
conditions 



To be concrete, let us consider the geometry shown in 
Fig. [1] and assume for the sake of simplicity that the 
strength of the SOI is zero in the leads. An incom- 
ing spinor valued wave is assumed to reach the device 
through lead /, and then it is generally partially reflected. 
Thus using expansion in terms of plane waves, we may 
write 



I*/} 



with E(k) = 



iE(k)i 



ikx 



/a(*) 



—ikx 



r 2 (k) 



dk, 



(24) 



Trlr 



= k 2 



(Note that e.g. the value of 



2m* h£l 

ka = 20.4 corresponds to a ring with a = 250 nm at the 
Fermi energy (11.13 meV) of InGaAs.) 

There is no incoming wave in the outgoing arm, that 

is 



/>oo 

I*//) = / e" 
Jo 



kx~iE(k)r ( tl(k) | ^ 
*2(fc), 



(25) 



In the interaction picture introduced by Eq. ([3]), the 
spinor components in the incoming and outgoing leads 



obtain time dependent phases, e.g, (fi(k), f 2 (k)) have to 
be replaced by (/i(fc)e i5T / 2 , f 2 {k)e~ li>T / 2 ). 

As a consequence of the presence of the leads, the index 
j of the eigenstates of the Hamiltonian |T]) does not need 
to be an integer any more, and in agreement with the 
notation used in Ref. (4lj we replace j by the continuous 
quantum number k. We expand the wave function in the 
two arms using this continuous variable and accordingly 
shift it from the index to the argument. In the interaction 
picture mentioned above we have 



l*«> 



B 



+ / B~{K)e 



u(K)e i ( K - 1 / 2 )*'' 



-iE~{k)t 



u(n)e 



i(K+l/2)<p 



dK. 
(26) 

A similar expression can be written for the spinor val- 
ued wave function in the lower arm, where the expansion 
coefficients will be denoted by C + (k) and C~(n). 

In these equations the functions fi(k) and f 2 (k) 
are assumed to be known, and we have to deter- 
mine ti(k),t2(k),ri(k), r 2 (k) and B ± {k),C ± {k) from the 
boundary conditions. For the sake of definiteness let us 
focus on junction 1. The requirement of continuity reads: 

\^!(x = 0,r)> = = 0,r)) = |* ; (^ 2 = 0,t)> , 

(27) 

and it has to be satisfied at all times r. Provided this 
holds, it can be shown that the requirement of vanishing 
energy current density at this junction is satisfied, if 



x=0 



dipi 



I* u (vi,t)> 



Vi=0 



dip 2 



(28) 



^2=0 



similarly to the case of a constant Hamiltonian, but now 
for arbitrary r. In other words, energy conservation fol- 
lows from the fact that the requirements ii) and iii) men- 
tioned in the introductory part of this section are satisfied 
at any given time. 

The set of equations (|2T[> and (|28|) together with 
the corresponding ones for junction 2 are solved most 
straightforwardly by Fourier transformation with respect 
to r. E.g, the first equation in l|27p is transformed as: 



|*/(2; = 0,r))e™ r dr 



\* u {i Pl = Q,T))j WT dr 



which leads to 



(29) 



/i(fc+M)\ (h{k+{w)) 

Mk-(w))J + \r 2 (k-(w)) 

B + (4(wm + (4(w))) + s+(4H)|v>+(4M)> 

B-(^(wj)\^(K^(w)))+B-(^(w)M-(^(w))) 



where k (w) — *Jw =p vjlja if the argument of the 
square root is not negative, and zero otherwise. Similarly, 
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Ki(w) and «2 i w ) are the two solutions of the equation 
E + {k) = w and similarly = £~(/t7/(u>)) = 



Additionally, 



B±(k) 



|^±(k)/9k| 



(30) 



and /, f also denote / and r divided by the modulus of 
dE{k)/dk. 

For a given value of w, Eqs. I[27p and (|28p and their 
counterparts at junction 2 leads to a closed set of 12 
linear equations with 12 unknowns (which are the ex- 
pansion coefficients r\,r2, etc. evaluated at certain val- 
ues of their respective arguments). This means that 
by sweeping w so that fc ± (w) covers the range where 
fi(k + (w)) and f2(k~(w)) are nonzero, we can obtain all 
nonzero values of the unknown functions. Thus inverse 
Fourier transformation can be used to calculate the solu- 
tion of the time dependent transport problem. (Note 
that the linearity of the equations implies that when 
fi(k + (w)) = f2(k~(w)) = 0, the solution also vanishes 
everywhere.) This method provides a general framework 
to investigate the conductance properties of the device 
with oscillating SOI strength. 



C. Discussion 

First we consider the case when the incoming wave has 
a narrow energy distribution, as shown by the dotted line 
in Fig. [5] For the sake of simplicity we assume a com- 
pletely unpolarizcd input spin state, and calculate the 
weights of the plane waves e zkx in the output. If the ex- 
ternal electric field were constant, this latter distribution 
of the spatial frequencies given by \ti(k)\ 2 + \t2(k)\ 2 would 
be zero outside the support of |/i(fc)| 2 + /2(A)! 2 , as an 
incoming plane wave with a given value of k would lead 
to transmitted states with the same wave number. (The 
linearity of the problem would forbid those values of k in 
the output, for which |/i(fc)| 2 + |/ 2 (fc)| 2 = 0.) However, 
as shown in Fig. [51 this is not the case, when the SOI 
strength oscillates. Besides the central peak correspond- 
ing to the wave numbers contained in the input, there 
are two additional values of k, where |ii(fc)| 2 + \t2(k)\ 2 
has pronounced maxima. Calculating the separation of 
the neighboring peaks in frequency units, we obtain that 
both are equal to the frequency of the SOI oscillations, 
v. Mathematically, this result is a consequence of the 
structure of the fitting equations, namely that e.g. fi(k) 
is connected to both t\{k) and t^yk 2 



via 2 ) 



h(k) 



h(^k 2 + u/a 2 ) 
t 2 (k) 

ti(fc) 

t 2 Wk 2 -u/a 2 ). 



(31) 



In a somewhat wider context, we can say that the effect 
is related to the time-energy uncertainty relation: when 
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22.50 
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FIG. 5: The weight |ti(fc)| 2 + \t 2 (k)\ 2 of the transmitted 
plane waves as a function of the wave number k (solid black 
curve). The incoming distribution |/i(fc)| 2 + |/2(fc)| 2 is shown 
by the dotted line. The tri plet struct ure is a consequence 
of scattering events k — > \Jk 2 ± D/a 2 . The parameters are 
A = 0.5, = 1.0. 



the characteristic time of a process is not infinitely long, 
it cannot correspond to a well-defined energy value. Let 
us also note that the "sideband currents" shown in Fig. [5] 
also appear when an oscillating scatterer is being placed 
in the rin g 54 i 55 . 

Considering the limit of infinitely narrow incoming dis- 
tribution, i.e, when fi(k) and f2(k) are either zero, or 
proportional to S(k — k Q ), we may ask what the conduc- 
tance of the device is. According to the previous results, 
in general, the transmitted (and reflected) state contains 
wave numbers ko and \/fc 2 ± vja 2 . (Note, however, that 
at zero temperature, if kg represents the Fermi wave num- 
ber, then the output corresponding to k = \fk^ — vja 2 is 
suppressed due to the occupation of levels with k < ko-) 
Then a direct calculation of the output current shows 
that generally it will not be constant, but oscillations 
with the frequency of the external field appear. Therefore 
we calculate average conductance, when time dependent 
cross terms disappear due to averaging over a period of 
T = 2tt/v. According to the Landauer-Buttiker formula, 
this averaged conductance in units of e 2 /% is given by 

G = ^|i 1 (fc )| 2 + |< 2 (A ; o)| 2 



yjk 2 + vja 2 



y/k 2 , - vja 2 
k 



h^kl + v/a 2 ) 



(32) 



t 2 Uk 2 -v/a 2 ) 



where the sum runs over two orthogonal inputs, 
e.g. (exp(ifcox), 0) and (0, e~xp(ikox)) . The result of this 
calculation is shown in Fig. [6] As we can see, the triplet 
structure shown in Fig. O appears again, here in the form 
of minima in the conductance. Close to the pronounced 
minimum of G around e.g. k m i n a = 20, there arc two 
additional minima corresponding to k = yfcj^jiti^?. 
This 'hole burning' in the transmission is the consequence 
of the fact that the minima of the third and fourth terms 
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FIG. 6: Time averaged conductance given by Eq. f!32 p for the 
case of A = 0.5 as a function of ka. 

in Eq. (|32j) are situated at different positions from those 
of the first two terms. In other words, scattering events 
that change the momentum, also modify the conductance 
of the device. 

Finally let us investigate the transmission of a wave 
packet through the ring. We choose here the case which 
in some sense is the opposite to the one discussed above, 
as now the distribution |/i(fc)| 2 + |/2(fc)| 2 is wide in k, 
resulting in a narrow wave packet in space. For the sake 
of defmiteness, the initial state is assumed to be spin- 
polarized in the positive z direction, and localized in the 
input arm with Gaussian envelope 

/ _ (*-*p)2 +ik \ 

l*OM = 0)>=l e ^ J (33) 

as shown in Fig.[7]for the case of a = 0.5a, xq — 1.5a and 
koa = 10. As we can see, when the wave packet reaches 
the ring, it is partially reflected, and the reflected waves 
interfere with the incoming packet, that leads to oscil- 
lating electron density in the input arm. The fraction 
of the wave packet that enters the ring travels along the 
two arms, interferes around the output junction and pro- 
duces an output wave packet. However, there is a certain 
probability for the electron not to leave the ring, there 
is a fraction of the wave packet that moves backwards 
toward the input junction, where interference phenom- 
ena can be observed again, and a weak reflected wave 
packet is formed that leaves the ring through the input 
junction. These processes are repeated periodically un- 
til the modulus of the wave function becomes negligible. 
As a consequence, both in transmission and reflection we 
can observe a series of wave packets with decreasing am- 
plitudes and increasing widths. The separation of these 
wave packets in time is essentially the round trip time of 
the wave packet in the ring, which is around ir/(akoQ). 
Effects related to the SOI are twofold: first the trans- 
mission probability (that determines the heights of the 
consecutive transmitted wave packets) depends on the 
amplitude and frequency of the SOI oscillations, and the 
spin directions in the ring are also determined by these 
parameters. 
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FIG. 7: Transmission of a packet through the quantum ring 
(A — 0.5, v = 0.5). The initial state is given by Eq. (|33[1 . with 
a = 0.5a, xo — 1.5a. The round trip time is given by r roun( j = 
7r/(fc a) = 0.2, where fc a = 10. See EPAPS Document No. 
[number will be inserted by publisher] for the related movie 
file. 

IV. RANDOM SCATTERERS, THERMAL 
FLUCTUATIONS AND AN APPLICATION 

In this section we investigate to what extent the ef- 
fects discussed so far are still present in a more real- 
istic context, when random scattering events and ther- 
mal fluctuations are also taken into account. To this 
end, first we introduce point-like random scatterers repre- 
sented by Dirac-delta potentials. That is, we add a term 
U(D) = u(D)S(4>) to the Hamiltonian {TJ), and solve the 
time dependent problem determined by this perturbed 
Hamiltonian. The strength of the potential, U(D), is 
random, it is drawn from a normal distributions, with 
zero mean and root-mean-square deviation D. The ran- 
dom fluctuation of the potential causes that the "aver- 
age" state of the system cannot be represented by a pure 
quantum mechanical state, it becomes a mixture that can 
be described by a density operator p{D). In practice, 
p(D) is calculated from several computational runs, with 
appropriate averaging. That is, for a particular value of 
the potential U(D), we obtain a solution spinor, that we 
write symbolically as |$[J7 (£))]) . When after M c compu- 
tational runs, the estimated density operator 

Pout(D) = ±. ]T MU(D)}) mU(D)}\ , (34) 

c n 

converges, we have all the possible information needed to 
describe what effects result from the disturbances char- 
acterized by the variable D. In this way, by tuning D we 
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can model weak disturbances (small D) as well as the 
case when frequent scattering events completely change 
the character of the transport process (corresponding to 
large values of D). 

The energy distribution of the input electrons can be 
taken into account by appropriate averaging over the pos- 
sible input energies. In thermal equilibrium at temper- 
ature T, the conductance of the device can be written 



G(T) = J p(E,T)G(E)dE. 



(35) 



^[1 + exp (E - E F )/k B T}- 1 and 



where p(E, T) = 
G(E) is given by Eq. (|32|) . (Note that this is essentially 
the Landauer-Buttikcr formula at finite temperature and 
low bias^.) For the sake of numerical convenience, we 
can convert the integral (|3"5"|) to a sum over the possible 
energies, meaning that the expression for p ut(T) is sim- 
ilar to Eq. (p4|) . but the weights of the projectors are not 
uniform, they are determined by the Fermi distribution. 

Considering first a closed ring, we found that an in- 
dividual scatterer shifts the energy levels according to 
the strength of the Dirac-delta potential. Consequently, 
when averaging over numerous random scattering events, 
we obtain a certain broadening of the possible energies as 
well as the Rabi frequencies. Investigating the time evo- 
lution of the expectation value S z , we found that for the 
parameters and initial state shown in Fig. [21 the height 
of the first maximum at t = 27rf2fl decreases when we 
increase D, the width of the random distribution. As 
a reference, we considered the case of time independent 
Hamiltonian (when ui(t) is constant in Eq. (Q])) and cal- 
culated the conductance Aharonov-Casher (AC) oscilla- 
tions for various values of D. We found that the visibility 



1(D) 



Gmax C^J G ra i n (^) 
Graax(D) + G m i n (D) 



(36) 



of these experimentally detectable oscillations decreases 
from unity to 1/2 while D increases from zero to D 1 / 2 = 
0.03 x Ep. (Ep denotes the Fermi energy.) At -Di/ 2 , the 
first maximum of S z is approximately 50% of the ideal 
value, thus a damped oscillation can be detected. On the 
other hand, the revival phenomena shown in Figs. [3] and 
5] are more sensitive to the broadening of the frequen- 
cies and consequently they are not expected to be de- 
tectable in usual samples. Note that this is in accordance 
with the expectations, as the "spintronic Schrodinger-cat 
sates" shown in Fig. 0] are highly nonclassical, and con- 
sequently they are exceptionally sensitive to any kind of 
environmental noised. 

The question to what extent the transport properties 
of the time dependent problem are modified by random 
scatterers and thermal fluctuations will be analyzed us- 
ing a possible application of our model. First we recall 
Figs. [5] and [5J showing a triplet structure in the trans- 
mission and the transmission probabilities for infinitely 



narrow incoming energy distribution. Clearly, the first 
effect is rather common when oscillating potentials are 
considered, but the second one is based on the specific 
spin dependent interference phenomena that appear in a 
ring structure. Combining these two effects, one can see 
that it is possible to find parameter values, where the ra- 
tio of the transmitted "direct" and "sideband" currents 
is strongly modified by the transmission profile. Par- 
ticularly, if the transmission probability for the central 
peak is suppressed by interferences (like around integer 
values of ka in Fig. [5]), sideband currents can provide 
stronger output signals than the direct one. (Usually, 
in other physical systems, it is mainly the strength of 
the oscillating potential that determines the relative in- 
tensity of the relevant frequencies in the output.) This 
situation is shown in Fig. [5] for the case of spin-up in- 
put in the z direction. Let us recall that while in the 
direct current this spin direction remains unchanged, the 
sideband current corresponds to oppositely polarized out- 
put spins. (Note additionally, that the positions of the 
minima of the transmission probabilities as a function of 
ka can be tuned by adding a constant SOI term to the 
Hamiltonian ([I}, a minimum can even be transformed 
into a maximum by using experimentally achievable SOI 
strengths 38 .) The question whether the two peaks are 
distinguishable, is clearly related to ratio of hi> and the 
width of the thermally broadened input energy distribu- 
tion: the shift induced by the oscillating SOI strength 
has to be larger than the width of the distribution. 
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FIG. 8: Direct (d) and sideband (s) currents when the input 
energy distribution (centered at Fermi energy corresponding 
to ka = 19.98) is broadened due to thermal fluctuations. The 
relevant part of the transmission profile seen for a wider in- 
terval in Fig. HJ] is also shown. In this case (parameters are: 
A — 0.75, v — 5.0) the sideband current can be considerably 
stronger than the direct one. The input spins are chosen to 
be polarized in the positive z direction and the transmission 
peaks has been normalized so that the maximum of the higher 
one is 1, for both temperatures, separately. 

Fig. [5] visualizes the case of low temperatures, and does 
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not take random scattering events into account. Now 
we investigate whether the effects is still visible when 
D ^ and the temperature is increased. Choosing ex- 
perimentally achievable temperatures^ Fig. [5] shows the 
ratio J I s (E)dE/ J Id{E)dE of the sideband and direct 
peaks (seen in Fig. |8f as a function of the strength D 
of the random scattering events. (Numerically, if Ep is 
assumed to be in the range of 10 meV, the temperatures 
corresponding to the curves shown in Fig. [9] have the or- 
der of magnitude of 100 mK.) As we can see, the effect of 
pronounced sideband peaks is still present for moderate 
values of D. 




Ct U i I , I i I i I , I i L 



0.000 0.002 0.004 0.006 0.008 0.010 0.012 

DIE 

F 

FIG. 9: The ratio of the integrated sideband and direct peaks 
shown in Fig. [8] as a function of the strength of the random 
scattering events. Integration was carried out with respect to 
the energy. 

Finally let us note that there are samples where the 
transport is due to many channels in the ring, thus our 
results are not directly applicable. For narrow rings, how- 
ever, the analysis of this section shows that some of the 
previously predicted effects - that are based on an ide- 
alized description - are stable enough against random 
scattcrcrs and thermal fluctuations to be visible. 



V. SUMMARY 

In this paper we investigated quantum rings with time 
dependent spin-orbit interaction (SOI). The effect of the 
sinusoidal SOI strength was first investigated in an iso- 
lated ring, i.e., a ring without any attached leads. In 



this case we have shown that for initial states with well- 
defined z component of the total angular momentum 
Rabi oscillations appear: The spin direction along the 
ring changes periodically in time. Considering the dy- 
namics of an initially localized wave packet, we demon- 
strated that appropriately chosen amplitudes of the SOI 
oscillations can lead to quasi-periodic time evolution: 
First the wave packet becomes delocalized (this process 
can be termed as collapse), then at the revival time it 
reassembles again. During this process (without thermal 
fluctuation or scattering events), there are time instants 
(partial revival times) when "spintronic Schrodinger-cat 
states" appear: superpositions of wave packets localized 
at different positions along the ring. 

In the second part of the paper we focused on the trans- 
mission properties of the device. We introduced a gen- 
eral method to treat the quantum mechanical scattering 
problem for a time dependent Hamiltonian. Based on 
this result, we have shown that in general the scattered 
outgoing plane wave states will have an energy distribu- 
tion which is different from that of the incoming states. 
For a single incoming plane wave with energy E, the out- 
put shall contain energy values E and E ± hv, that is, 
shifts corresponding to the frequency of the SOI oscil- 
lation appear. If the input is spinpolarized e.g. in the 
positive z direction, the sideband current corresponding 
to E + hv is related to oppositely polarized output. Ad- 
ditionally, it was shown that one can find parameter val- 
ues, where sideband currents can provide stronger out- 
put signal than the direct one. This transport property 
is specific to quantum rings: sideband currents appear 
in a wide class of physical systems, but the most impor- 
tant role in the suppression of the direct current is played 
by spin dependent quantum interference, which is closely 
related to the geometry of the device. Additionally, this 
effect is still visible at finite temperatures even if random 
scattering events modify the the dynamics. According 
to our calculations, Rabi oscillations can also survive a 
moderate level of temperature induced fluctuations. 
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